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The whispering gallery modes (WGMs) of quartz microspheres are investigated for the purpose 
of strong coupling between single photons and atoms in cavity quantum electrodynamics (cavity 
QED). Within our current understanding of the loss mechanisms of the WGMs, the saturation 
photon number, no, and critical atom number, No, cannot be minimized simultaneously, so that 
an "optimal" sphere size is taken to be the radius for which the geometric mean, \/no x No, is 
minimized. While a general treatment is given for the dimensionless parameters used to characterize 
the atom-cavity system, detailed consideration is given to the D2 transition in atomic Cesium at 
Ao = 852nm using fused-silica microspheres, for which the maximum coupling coefficient g a /(2-K) m 
750MHz occurs for a sphere radius a — 3.63um corresponding to the minimum for no ~ 6.06 x 10~ 6 . 
By contrast, the minimum for No ~ 9.00 x 10 -6 occurs for a sphere radius of a — 8.12/im, while the 
optimal sphere size for which \/no x No is minimized occurs at a = 7.83/xm. On an experimental 
front, we have fabricated fused-silica microspheres with radii a ~ 10pm and consistently observed 
quality factors Q > 0.8 x 10 7 . These results for the WGMs are compared with corresponding 
parameters achieved in Fabry- Perot cavities to demonstrate the significant potential of microspheres 
as a tool for cavity QED with strong coupling. 

PACS numbers: 42.55.Sa, 42.50.Ct, 32.80.-t 
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I. INTRODUCTION 

Motivated by the pioneering work of Braginsky and 
Ilchenko (!]], some of the highest quality optical res- 
onators to date have been achieved with the whis- 
pering gallery modes (WGMs) of quartz microspheres 
||] . Over the wavelength range 630 — 850nm, qual- 
ity factors Q w 8 x 10 9 have been realized, and cav- 
ity finesse T = 2.3 x 10 6 demonstrated [§, [§. Such 
high quality factors make the WGMs of small dielec- 
tric spheres a natural candidate for use in cavity QED 
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While much of the work regarding quartz microspheres 
has centered around achieving the ultimate quality fac- 
tors D , the quality factor of the resonator is but one of 
the factors that determines the suitability of the WGMs 
for investigations of cavity quantum electrodynamics in a 
regime of strong coupling. In this case, the coherent cou- 
pling coefficient, g, for a single atom interacting with the 
cavity mode must be much larger than all other dissipa- 
tive rates, including the cavity decay rate, k, and the rate 
of atomic spontaneous emission, 7; namely g (k, 7). 
Note that 2g — Q, gives the Rabi frequency associated 
with a single quantum of excitation shared by the atom- 
cavity system [Ml |2^]. The atom- field interaction can 

be characterized by two important dimensionless param- 

2 

eters: the saturation photon number, no a ^, and the 
critical atom number, Nq oc ^f. Since these parame- 
ters correspond respectively to the number of photons re- 
quired to saturate an intracavity atom and the number of 
atoms required to have an appreciable effect on the intra- 
cavity field, strong coupling requires that (no, A^o) <C 1. 
Ideally one would hope to minimize both of these param- 
eters in any particular resonator. Unfortunately, within 



the context of our current understanding of the loss mech- 
anisms of the WGMs ||, the critical parameters (n 0l N ) 
cannot be minimized simultaneously in a microsphere. 

Motivated by these considerations, in this paper we ex- 
plore possible limits for the critical parameters (no,iVo) 
for the WGMs of quartz microspheres. Following the 
analysis of Refs. J|, ||, |l8) , we study the particular case of 
a single atom coupled to the external field of a WGM near 
the sphere's surface. We show that there are radii that 
minimize (no, No) individually, and that there is an "op- 
timal" sphere size that minimizes the geometric mean, 
V^o x No, of these two cavity-QED parameters and al- 
lows both parameters to be near their respective minima. 
We also report our progress in the fabrication of small 
microspheres with radii a ~ 10/im, and compare our ex- 
perimental results for Q with those from our theoretical 
analysis. Finally, we present a detailed comparison for 
the state of the art and future prospects for achieving 
strong coupling in cavity QED for both microsphere and 
Fabry-Perot cavities. Throughout the presentation, we 
attempt to develop a general formalism that can be ap- 
plied to diverse systems. However, for definitencss we 
also present results for a particular system of some inter- 
est, namely an individual Cesium atom coupled to the 
WGMs of quartz microspheres. 



II. MODES OF A MICROSPHERE 

Solving for the mode structure of the resonances of a 
dielectric sphere in vacuum is a classic problem in elec- 
tricity and magnetism, and the resulting field distribu- 
tions have been known for some time [{23| . The electric 
field of the TM, electric type, modes inside and outside a 
sphere of refractive index n at free-space wavelength Aq 
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are respectively, 
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where a is the radius of the sphere, k — is the wave 
vector inside the sphere, ji (x) is the spherical Besscl func- 
tion, hi(x) is the spherical Hankel function, (r, 9, <j>) are 
unit vectors, and the ' refers to differentiation with re- 
spect to the argument. Note that the TM modes have a 
predominantly radial electric field vector. 

In order to satisfy the boundary conditions at the sur- 
face of the microsphere, the tangential components of the 
mode function immediately inside and outside the sphere 
must be equal. However, there is a discontinuity in the 
radial component of the electric field at the dielectric 
boundary (as can be seen from Fig. [j].) The eigenmodes 
are determined by solving for the roots of a characteristic 
equation |23| , which can be reduced to 
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Throughout this paper, we normalize the mode func- 
tions such that their maximum value is unity. This con- 
dition then yields for the I = m modes of the sphere 

#in(r,M)= N(l + l) J -^-sm l (9)e^r 
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FIG. 1: The magnitude of the normalized mode function as a 
function of radius for the TM mode of a 5.305/im radius sphere 
(p = 1, 1 — m — 50) with 6 — ^ and <j) — for a wavelength of 
Ao = 852.359nm and index of refraction n — 1.45246. In our 
case, the function is normalized to have a maximum value of 
unity. Note that there is a discontinuity at the surface. 
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and N is the normalization factor. Because we will re- 
quire the field outside the sphere to be as large as possi- 
ble, we will choose the p = 1 modes. Also, because the 
coherent coupling constant g oc 



where Vp is the 



cavity mode volume, we choose the / = m modes, since 
they yield the smallest electromagnetic mode volume, as 
will be explained in the next section. 



III. ELECTROMAGNETIC MODE VOLUME 



The effective mode volume Vp associated with the elec- 



tromagnetic field distribution ^(r, 9, c 



is given by 



dV 



where 



e(r) = 



n if r < a, 
1 if r > a. 



(9) 
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and P corresponds to the (p, I, m) mode. Vq is the quan- 
tization volume discussed in Ref . [Q . As long as a radius 
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tq is chosen large enough to include the effects of the 
evanescent field, the mode volume is relatively insensi- 
tive to the particular choice of quantization radius [^3. 
As discussed more extensively in Refs. [ 
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the 



teraction between the internal atomic degrees of freedom 
and the intracavity field is characterized by the coherent 
coupling constant g(r, 9, </>), where 



g(r,d, ( f>) = g o ^ l ' m \r,9,0) 



(11) 



and 



Fig. ||) . One might at first believe that this value for the 
radius represents the optimal sphere size for use as a cav- 
ity with single atoms. However, while the mode volume 
Vp plays an important role in determining the coupling 



constant (Eq. 12), it is not the only parameter relevant 



to cavity QED with single atoms in a regime of strong 
coupling. As discussed in the next sections, the qual- 
ity factor, Q, of a WGM has a strong dependence on the 
sphere radius, and must also be considered in an attempt 
to optimize the critical atom and saturation photon num- 
bers. 
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Note that in the absence of damping, 2g (r) gives the fre- 
quency for Rabi nutation associated with a single photon 
in the cavity for an atom initially in the ground state lo- 
cated at position r within the mode. Therefore, in order 
to maximize the coupling strength, one must endeavor to 
minimize the cavity mode volume. 

In order to derive an answer that can be applied to dif- 
ferent wavelengths, one can define a dimensionless mode 
volume parameter, V, and plot as a function of a dimen- 
sionless sphere size parameter, x, defined as: 
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where Vp is the cavity mode volume, n is the index of 
refraction at the free-space wavelength Ao, and a is the 
sphere radius. The plots then only depend on the index 
of refraction (see Fig. ||). 

Naively, one might assume that the sphere should be 
made as small as possible in order to minimize the elec- 
tromagnetic mode volume, and hence to provide a max- 
imum for go and hence globally for g(r). However, as 
shown in Figs. || and |[ the mode volume for the TM 
modes of a quartz microsphere actually passes through 
a minimum at some particular radius ao. This behav- 
ior can be understood by noting that for a < ao, the 
intrinsic, radiative losses are increasing rapidly and ul- 
timately cause the mode to no longer be well-confined 
by the sphere, with a concomitant increase of the mode 
volume. Note that in Fig. || and subsequent figures, we 
give results for n ~ 1.45 corresponding to fused silica, 
as well as for n = 2.00 and n — 3.00. These latter 
cases serve to illuminate the role of n as well as being 
applicable to other materials (i.e., the index of refrac- 
tion for GaAs is n — 3.4 for A = 1550nm (25)). For 
a very low-OH fused silica microsphere at Ao = 852nm 
(the wavelength of the Di transition in atomic Cesium) 
with index of refraction n — 1.45246, the minimum mode 
volume Vp un s» 28.4/im 3 occurs for radius a as 3.73/im 
corresponding to mode numbers p — 1 , 1 = m = 34 (see 



IV. LOSSES IN DIELECTRIC SPHERES 

For fused silica spheres with radius a > 15//m, the 
effect of intrinsic radiative losses can be safely neglected, 
since they allow quality factor Q > 10 21 , as illustrated 
in Fig. ||. Such large values of Q greatly exceed those 
imposed by technical constraints of material properties, 
such as bulk absorption and surface scattering. 

However, as one moves to very small spheres with ra- 
dius a < 10/xm, the intrinsic radiative Q falls steeply 
enough to become the dominant loss mechanism even in 
the face of other technical imperfections. When assessing 
the usefulness of microspheres for cavity QED, one must 
account for the entire set of loss mechanisms to deter- 
mine the optimal size for the microsphere, which is the 
subject to which we now turn our attention. 

The quality factors of the WGMs of fused silica micro- 
spheres arc determined by several different loss mecha- 
nisms. The overall quality factor can then be calculated 
by adding the different contributions in the following way 



Q' 1 = Q^ + Q 
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mat 3 
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Qmat — Qs.s. + Qw + Qbu) 



lk- 



(15) 
(16) 



where Q la d is due to purely radiative losses for an ideal 
dielectric sphere and Q m at results from non-ideal mate- 
rial properties. The principal mechanisms contributing 
to Q m at are scattering losses from residual surface in- 
homogeneities (Q s .s.), absorption losses due to water on 
the surface of the sphere (Q w ), and bulk absorption in 
the fused silica (Qbuik)- The intrinsic material losses are 
known very accurately, since they arise from absorption 
in the material at the wavelength of concern ^(|. Con- 
siderably greater uncertainty is associated with the losses 
due to surface scattering and absorption due to adsorbed 
material on the surface of the sphere, of which water is 
likely the principal component. Wc will adopt the mod- 
els for these losses presented in Refs. ^ ||, extrapolated 
to the regime of small spheres of interest here. 



A. Intrinsic Radiative Losses 

The contribution to the quality factor for purely ra- 
diative effects, Q ra d, can be derived by following the ar- 
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FIG. 2: fa) The dimensionless volume parameter, V (defined 
by Eq. h3j), as a function of the dimensionless size parameter, 
x (defined by Eq. p"l| ). The solid line is for an index of re- 
fraction n = 1.45246, the index of refraction for fused silica at 
A = 852nm, with a minimum of V = 34883.4 for x = 39.9469 
(I — m — 34). The dotted line is for an index of refraction 
n = 2.00, with a minimum of V = 15596.2 for x = 18.9864 
(I — m — 14). The dashed line is for an index of refraction 
n = 3.00, with a minimum of V = 11546.4 for x = 10.2748 
(I = m = 6). (b) Because the index of refraction for fused 
silica varies from n = 1.444 at Aq = 1550nm to n = 1.458 
for Ac, = 600nm (see Fig. §), this plot of the dimensionless 
volume parameter, V, as a function of the dimensionless size 
parameter, x, is made for that range of values. The solid line 
is for an index of refraction n = 1.44, with a minimum of 

V = 36247.5 for x = 40.9812, (I = m = 35). The dotted 
line is for an index of refraction n = 1.45, with a minimum 
of V = 35161.1 for x = 41.0036, (Z = m = 35). The dashed 
line is for an index of refraction n — 1.46, with a minimum of 

V = 34129.1 for x = 39.9631, (l = m = 34). 
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FIG. 3: The electromagnetic mode volume, Vp, for the TM 
modes of a very low-OH fused silica microsphere as a function 
of sphere radius at the wavelength Ao = 852nm for the D2 
line of atomic Cesium. The minimum, 28.4/im 3 , occurs for 
radius ao ~ 3.73/im corresponding to mode numbers p = 1 
and I = m = 34. 
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FIG. 4: The index of refraction of very low-OH fused silica as 
a function of wavelength. 



Lorenz-Mie theory @. We find from Ref. [g7| that 

Qr ad = Ul + l)n^ b (n 2 -l) 1/2 e^, (17) 



where 



Ti = ( I + - ) (rji - tanh?7;) , 



(18) 



guments presented in Ref. |27 . These losses are due to 
the leakage of light from the resonator due to its finite 
dielectric constant and radius of curvature. The results 
can then be compared to numerical results obtained by 
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FIG. 5: (a) Semi-log plot of the radiative quality factor, Q ra d, 
for various indices of refraction as a function of the dimen- 
sionless size parameter, x = 2 T n " . (b) Semi-log plot of the 
radiative quality factor, Q ra d, as a function of sphere radius 
for a wavelength of Ao = 852.359nm (index of refraction is 
n = 1.45246). 




reaches 1% for I — 18. However, the error is less than 
0.2% for I = 76 (This is the optimal sphere size discussed 
in Section |Vl| ) . Fortunately, the expressions are valid in 
the regimes for which we are concerned. This has been 
confirmed by making comparisons with numerical values 
obtained using Lorenz-Mie scattering theory. 

From Fig. a, we see that the radiative Q falls approxi- 
mately exponentially as the radius a is decreased, and can 
become quite important as the sphere size is decreased 
below 10/im. For example, for a 15/im radius sphere and 
a wavelength A = 852.359nm, Q rad «2x 10 21 . There- 
fore, the net quality factor would most certainly be dom- 
inated by other loss mechanisms in Eq. [l^. However, for 
a 7jUm radius sphere, Qrad ~ 4 x 10 s , and the radiative 
losses can play a crucial role in the characteristics of the 
spheres that are optimal for use in cavity QED. 



B. Material Loss Mechanisms 

The quality factor due to bulk absorption, Qbuik, in 
fused silica is actually known very well, since this depends 
only on the absorption of the material at the wavelength 
of concern 101: 



_ 27m 
Wbulk — — r~J 



(22) 



where n is the index of refraction, and a is the ab- 
sorption coefficient of the material. From Fig. ^ we see 
that for very low-OH fused silica, the absorption coeffi- 
cient at 852nm is a w 4.5 x 10~ 4 m _1 Q. This would 
correspond to a quality factor of Qbuik ~ 2.4 x 10 10 . 
Fused silica has a minimum in its absorption coefficient 



of aw 1.5 x 10" 
factor of Qbuik r 



J m at 1550nm, which yields a quality 
3.8 x 10 



li 



The quality factor due to surface scattering, Q s .s., and 
absorption by adsorbed water, Q w , has also been studied 
and modelled, albeit for larger spheres with a > 600/im. 
For losses due to surface scattering, we follow the work 
of Refs. 10. 01 and take 



and 



1 + 



TE modes, 

1 TM modes. 



(20) 



(21) 



Also, n is the index of refraction and t° is the p th zero 
of the Airy function Ai. This p corresponds to the mode 
number (p,l,m). In our case, we are only interested in 
the p = 1 modes of the sphere to maximize the elec- 
tromagnetic field outside the sphere while maintaining a 
small mode volume. Note that these expressions for Q ra d 
become invalid in the limit of small / mode numbers. The 
error in the mode functions used to derive these results 
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where e = n 2 is the dielectric constant and aB ~ 5nm 2 
is an empirical parameter determined by the size and 
correlation length of the distribution of residual surface 
inhomogeneities. This quantity was reported in Ref. ^ 
based upon atomic force microscopy measurements of a 
microsphere. 

The quality factor due to water adsorbed on the sur- 
face, Q w , is given by H 
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where 5 ~ 0.2nm is an estimated thickness for the water 
layer, and j3 w ~ 4.33m _1 is the absorption coefficient of 
water at 852nm. 
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FIG. 6: The quality factor, Qbuik, from Eq. |22j for a very 
low-OH fused silica microsphere as a function of wavelength. 
Because fused silica has a minimum in absorption at 1550nm, 
there is a maximum for the quality factor due to bulk absorp- 
tion of Qbuik ~ 3.8 x 10 11 . At 852nm, the quality factor due 
to bulk absorption is Qbuik ~ 2.4 x 10 10 . 

Combining these various results, we display in Fig. |?] 
a curve for the quantity Q ma ± as a function of sphere 
radius, a, for a wavelength Ao = 852nm. This same fig- 
ure shows the quality factor, Q ra( j, set by intrinsic ra- 
diative losses (Eq. [Tl]), as well as the overall quality 
factor, Q — Jp rad , < ^? lat . From this plot, we see that the 
radiative losses dominate the overall quality factor below 
a radius of a < 8/im, while the losses due to material 
properties are most significant for a > 8/im. Because of 
the extremely steep dependence of <5 ra d on sphere size, 
the point of transition from material to radiative domi- 
nated loss should be reasonably insensitive to details of 
the models employed to describe the material losses. Al- 
though we focus our attention here on the wavelength 
appropriate to the particular case of the Di transition in 
atomic Cesium, a similar analysis could be carried out 
for other wavelengths of interest using the above formal- 
ism, as for example the 25 — * IP transition at 1.083/im 
in metastable Helium. 



V. THE STRONG COUPLING REGIME 
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FIG. 7: (a) Semi-log plot of the quality factors due to the 
various loss mechanisms discussed in Section |^j for a very 
low-OH fused silica microsphere as a function of sphere radius 
for the I — m, TM modes at a wavelength of Ao = 852nm. 
In particular, traces are shown for the quality factor due to 
purely radiative losses (Q ra d), the bulk absorption of fused 
silica (Qbuik), both radiative losses and bulk absorption, the 
three loss mechanisms comprising Q mat : (Qbuik, Qs.a., Qw), 
and the predicted Q due to all four loss mechanisms, (b) 
This linear plot zooms in on the region of interest at the 
transition where the radiative losses become the dominant 
loss mechanism. The plot contains the quality factor due 
to purely radiative losses (Q ra d), the three loss mechanisms 
comprising Q mat : (Qbuik, Qs.s., Qw), and the predicted Q due 
to all four loss mechanisms. 



The ultimate goal that we consider here is to employ 
the WGMs of quartz microspheres as cavity modes for 
achieving strong coupling to atoms within the setting of 
cavity QED. The atom of choice in this paper is Cesium, 
and in particular, the D2 (F = 4 ^ F' — 5) transition 
in Cesium at Ao = 852.359nm as an illustrative exam- 
ple. Such an analysis allows a direct comparison with 
the state of the art in Fabry-Perot cavities [^9| . 

The coupling coefficient g(f) is the coupling frequency 
of a single atom to a particular cavity mode and cor- 
responds to one-half the single photon Rabi frequency 



pl[ p2[ . For an atom located just at the outer surface 
of the microsphere (i.e., in vacuum) and interacting with 
a whispering gallery mode P — (p,l,m), the coupling 
coefficient is given by B 



g(a) = g a = ~/± ^out(a) 




(25) 



where a is the sphere radius, 3p- 



2.61MHz is the trans- 



verse spontaneous decay rate for our transition in Ce- 
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sium, Vq 



3cAg 

47T7 J_ 



is the effective volume of the atom for and 



purely radiative interactions, and Vp is the electromag- 
netic mode volume of the whispering gallery mode des- 
ignated by P — (p, I, m). 

Armed with a knowledge of g, we are now able to de- 
termine certain dimensionless parameters relevant to the 
strong coupling regime of cavity QED. In particular, we 
consider an atom-cavity system to be in the strong cou- 
pling regime when the single-photon Rabi frequency, 2g, 
for a single intracavity atom dominates the cavity field 
decay rate, k, the atomic dipole decay rate, and the 
inverse atomic transit time, T^ 1 |2l], ^2|. We will defer 
further discussion of T , however, this requirement re- 
lates to the need for atomic localization [Q, In the 
strong coupling regime, important parameters for char- 
acterizing the atom-cavity system are the two dimension- 
less parameters: the saturation photon number, no, and 
the critical atom number, Nq. The saturation photon 
number, given by 



n 



2.9 2 



(26) 



corresponds to the number of photons required to sat- 
urate an intracavity atom p2|| . The critical atom 
number, defined by 



No 



2t_la 

g 2 



(27) 



corresponds to the number of atoms required to have an 
appreciable effect on the intracavity field |2^] . Ideally, 
one hopes to minimize simultaneously both the critical 
atom number, Ao, and the saturation photon number, 
no, which corresponds to simultaneous maxima for both 
-2— and 4- 

The saturation photon number and critical atom num- 
ber are useful because of their physical meaning. How- 
ever, one can define a new dimensionless parameter 



p= ***V P 1 



3Ag 



*out(a) 



(28) 



that corresponds to the cavity mode volume in units of 
A 3 weighted by the inverse of the strength of the mode 
function at the atomic position. This enables the equa- 
tions for the saturation photon number and critical atom 
number to be expressed as: 



and 



where 



n 



Nq 



p 



4 Q atom 

P 

Q cavity 
TTC 



A071 



(29) 



(30) 



(31) 



n ■ - —— 

^cavity — 



TTC 

X k' 



(32) 



This parameter, [3, then also determines the coupling co- 
efficient in the following manner: 



9W 



2-7TC71 

~p\ 



(33) 



Therefore, we see that one can use a single parameter, (3, 
combined with the properties of the atom to be used (Ao 
and 7_i_) and the quality factor of the resonator, Qcavity, 
to determine the three parameters (no, Nq, g ) of impor- 
tance in determining the quality of an atom-cavity sys- 
tem. 

Figs. H and ^| are plots of this dimensionless parameter 
P and of -js as functions of the dimensionless size pa- 
rameter x = 2 T" a for a few values of index of refraction. 

An 

Because the index of refraction for fused silica varies from 
n = 1.444 at A = 1550nm to n = 1.458 for A = 600nm 
(see Fig. ||), Figs. 0b andj§b are made for that range of 
values. From Figs.^ and Bone sees that there is a mini- 
mum for (3 and a maximum for that depends on the 
index of refraction. 



VI. STRONG COUPLING WITH CESIUM 

The results of the previous section can now be used to 
determine the saturation photon number, no, the criti- 
cal atom number, Ao, and the coupling coefficient, 17(a), 
for any atomic transition. In our case, we are concerned 
with the D2 transition in Cesium (Ao = 852.359nm). For 
this transition, the spontaneous transverse decay rate is 
2.61MHz. Also, at this wavelength the index of re- 



2- 



fraction for fused silica is n = 1.45246. 



This allows one to 



compute the coupling coefficient, g(a) — y ■ Fig. 

shows that there is a maximum of #- = 749.986MHz 
for a radius a — 3.63/im, (/ = m = 33). Interestingly, 
because we are restricted to having the atom couple to 
the external field of the microsphere, the maximum in 
the coupling coefficient, 5(a), does not coincide with the 
minimum for the mode volume, Vp (see Figs. || and fl0| .) 

The saturation photon number, no, is proportional to 
the dimensionless parameter as shown in Eq. 



29. Since 



the factor of proportionality is a constant that depends 
only on the properties of the particular atom of concern, 
the curve is determined by that of (3 along with the qual- 
ity factor of the atomic resonance (in our case Cesium), 
which is given by Eq. [3l] to be Q a tom = 6.738 x 10 7 . 
Fig. [ll] is a plot of the saturation photon number for 
the D2 transition in Cesium as a function of sphere size. 
Fig. [ll] shows that there is a minimum for the satura- 
tion photon number of no = 6.05527 x 1CP 6 for a sphere 
radius of a = 3.63163/xm (I = m = 33). 

The critical atom number, Nq, is also proportional to 
the dimensionless parameter (3 as shown in Eq. 30 . How- 



ever, its factor of proportionality is the quality factor of 
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FIG. 8: (a) The dimensionless parameter /3 as a function of 
the dimensionless size parameter x = 2 ™° . For an index of 
refraction n = 1.45246 (i.e., the index of refraction for fused 
silica at Ao = 852nm), there is a minimum of /3 = 1632.01 
for x = 38.8833, (I — m = 33). For an index of refraction 
n = 2.00, there is a minimum of f) = 221.124 for x = 17.8763, 
(I — m = 13). For an index of refraction n = 3.00, there is 
a minimum of /3 = 45.3744 for x = 10.2748, (I = m = 6). 
(b) Because the index of refraction for fused silica varies from 
n = 1.444 at A = 1550nm to n = 1.458 for Ao = 600nm (see 
Fig. ^]), this plot is made for that range of values. For an index 
of refraction n = 1.44, there is a minimum of (3 = 1753.92 
for x = 39.9188, (I — m — 34). For an index of refraction 
n = 1.45, there is a minimum of /3 = 1653.7 for x — 38.8778, 
(I — m — 33). For an index of refraction n = 1.46, there is a 
minimum of = 1561.45 for x = 37.8348, (I = m = 32). 



FIG. 9: (a) The dimensionless parameter -4g as a function of 
the dimensionless size parameter x = 2 ^" a . For an index of 
refraction n — 1.45246 (i.e., the index of refraction for fused 
silica at Ao = 852nm), there is a maximum of = 0.0247536 
for x, — 38.8833, (I = m = 33). For an index of refraction n — 
2.00, there is a maximum of ^ = 0.0672484 for x = 17.8763, 
(I — m = 13). For an index of refraction n = 3.00, there is 
a maximum of ^ = 0.148455 for x = 10.2748, (I = m = 6). 
(b) Because the index of refraction for fused silica varies from 
n = 1.444 at Ao = 1550nm to n = 1.458 for Ao = 600nm (see 
Fig. 0), this plot is made for that range of values. For an index 
of refraction n = 1.44, there is a maximum of = 0.0238779 
for x — 39.9188, (I — m = 34). For an index of refraction n — 
1.45, there is a minimum of = 0.0245908 for £ = 38.8778, 
(I = m — 33). For an index of refraction n = 1.46, there is a 
minimum of ^ = 0.0253068 for x = 37.8348, (I = m = 32). 



the resonator, Q 'cavity, which has a very strong depen- 
dence on the sphere radius, a, in the region below 10/im 
(see Fig. Q). Therefore, the minimum for the critical 
atom number does not occur for the same sphere size 
as for the saturation photon number. Fig. |l2| is a plot 
of the critical atom number as a function of sphere size. 
Using for Q ca vity the model that incorporates all of the 
loss mechanisms discussed in section |V] (radiative losses, 



bulk absorption, surface scattering, and absorption due 
to water on the surface), we find that the minimum for 
the critical atom number No — 8.99935 x 10 -6 occurs for 
a sphere radius of a — 8.12015/im (I — m — 79). At this 
radius, the coupling coefficient is #- = 304.16MHz. 

Unfortunately, as illustrated in Fig. [l3|, the minima 
for the two parameters, no and No, do not occur for the 
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FIG. 10: The coupling coefficient, as a function of sphere 
size for the D2 transition in Cesium (Ao = 852.359nm). There 
is a maximum of ^- = 749.986MHz for a sphere radius of 
a = 3.63163/im, (I = m = 33). Note that the maximum for 
does not coincide _with the minimum for the cavity mode 



volume, Vp (see Fig. 
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FIG. 11: The saturation photon number, no, as a func- 
tion of sphere size for the D2 transition in Cesium (Ao = 
852.359MHz). There is a minimum no = 6.05527 x 10" 6 for a 
sphere radius of a = 3.63163/xm (I = m = 33). At this radius, 
the coupling coefficient is = 749.986MHz. 



same sphere radius. However, if one uses the minimum 
of the geometric mean of the two parameters, each can 
have a value near its respective minimum. The mini- 
mum of the geometric mean occurs for a sphere radius 
a = 7.83038/im (I = m = 76). For this sphere size, the 
coupling coefficient is j- = 318.333MHz, the saturation 
photon number is no = 3.36107 x 10 -5 , and the critical 
atom number is N = 9.27834 x 10~ 6 . Therefore, each 
cavity QED parameter can be made to achieve simulta- 
neously a value near its respective minimum. 
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FIG. 12: The critical atom number, Ao, as a function 
of sphere size for the D2 transition in Cesium (Ao = 
852.359MHz). There is a minimum iV = 8.99935 x 10" 6 
for a sphere radius of a = 8.12015/im (I = m = 79). At this 
radius, the coupling coefficient is = 304.16MHz. This plot 
of the critical atom number incorporates the model for the 
quality factor of the resonator, Q cavity, outlined in section 
tVl for the four loss mechanisms: bulk absorption, surface 
scattering, absorption due to water on the surface, and radia- 
tive losses. The dark grey region is bounded by the effects of 
purely radiative losses. The light grey region is bounded by 
the effects of both radiative losses and bulk absorption. 



VII. PROGRESS IN SMALL SPHERE 
MANUFACTURE 



A large portion of the work being done on microspheres 
has been to push the quality factors of the spheres to 
record levels [p], || . This effort has produced some of the 
highest finesse (J- = 2.3 x 10 6 ) optical cavities to date 
with quality factors Q ~ 10 10 ||, j|. However, we have 
seen that Q is not the only relevant factor in determining 
the suitability of the WGMs for cavity QED in a regime 
of strong coupling. In general, the preceding analysis 
demonstrates the requirement to push to microspheres of 
small radius, a < 10/im. Unfortunately, the experiments 
that have achieved the highest quality factors and which 
have investigated certain material loss mechanisms are 
of rather larger size, and hence not optimal for cavity 
QED in a regime of strong coupling. For example, the 
experiment of Ref. § achieved a quality factor of Q = 
7.2 x 10 9 at 850nm in a sphere of radius a = 340/im. 

To explore the possibilities of cavity QED with strong 
coupling in substantially smaller spheres, we have un- 
dertaken a program to study fabrication techniques for 
quartz microspheres with a < 30/im, while still maintain- 
ing high quality factors. We have been able to fabricate 
10/im radius spheres using an oxygen- hydrogen micro- 
torch to melt the ends of very low-OH fused silica rods 
to form a sphere on the end of a stem. Light is then 
coupled to the sphere using frustrated total internal re- 
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FIG. 13: This plot shows the two parameters, (no, No), of 
importance to cavity QED as a function of sphere radius. 
The geometric mean of these two parameters is also plotted. 
The solid line represents the saturation photon number, no, 
the dashed line gives the critical atom number, No, and the 
dotted line shows the geometric mean of the two parameters, 
\/no x No. The minimum of each plot corresponds to the 
following dimensionless parameters: no = 6.05527 x 1CF 6 for 
a = 3.63163Atm (I = m = 33), and iV = 8.99935 x 1(T 6 
at a = 8.12015/mi (I — m — 79). The two curves cross at 
a = 7.03/im with n = No = 2. 56 x 1CT 5 . The geometric 
mean of these two parameters, v'no x No, is minimized for 
a — 7.83038um (I = m = 76) . For this radius, the parameters 
are: n = 3.36107 x 10~ 5 and iV = 9.27834 x 10" 6 . Note that 
the curve for No assumes the model for the Q discussed in this 
paper, and that the coupling coefficient g (f) is evaluated at 
the maximum of the mode function for r — a. 



flection of a prism, as in Refs. ||, |[ |3lJ . Our observations 
demonstrate that spheres of this size can be made consis- 
tently to have quality factors Q > 0.8 x 10 7 . While this 
is encouraging progress, the resulting Q is two orders of 
magnitude smaller than the theoretical maximum of ap- 
proximately 1.3 x 10 9 for this size based upon the model 
discussed in Section IV. 



One possible reason for this discrepancy could be the 
importance of minimizing the ellipticity of the small 
spheres. Because the small resonators fabricated by our 
technique have a stem protruding out of them, they are 
far from spherical. When coupling to an I = m mode in 
spheres with a > 100/im and hence large I, the mode is 
tightly confined to the equator; therefore, the poles do 
not have an appreciable impact on the mode structure or 
quality factor. In this case, it is not of critical importance 
to have the best sphere possible, but rather the best great 
circle possible to achieve large quality factors. However, 
this is not the case in small spheres with a < 10/im. As a 
decreases, the I — m modes occupy an increasingly larger 
proportion of the sphere in polar angle, and the elliptic- 
ity of the sphere becomes increasingly important in de- 
termining the mode structure as well as the Q. However, 
while there is certainly room for improvement in our fab- 
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FIG. 14: The solid line gives a parametric plot of the criti- 
cal atom number, iVo, and the saturation photon number, no, 
for fused-silica microspheres and the D2 transition of atomic 
Cesium, incorporating the loss mechanisms outlined in sec- 
tion |vl The dark grey region is bounded by the effects of 
radiative losses. The light grey region is bounded by the ef- 
fects of bulk absorption and radiative losses. This plot also 
offers a comparison of experimental and theoretical cavity 
QED parameters for microsphere and Fabry-Perot cavities. 
■ represents the current state of the art for cavity QED in 
Fabry-Perot cavities as in Ref. po) |. □ is a projection of the 
practical limit for Fabry-Perot cavities based upon Ref. fjjl . 
▲ represents the 60/im radius sphere implemented for cavity 
QED in Ref. [^lj. T is the current state of the art in 10/im 
microspheres based upon the results presented in section [VlL 
A is the currently achievable Q with the optimal sphere size 
of 7.83/im based upon the analysis of sections IV and ^ V 
is the theoretically achievable Q ~ 9.76 x 10 s at the optimal 
sphere size, a w 7.83/xm. 



rication technique and in the resulting mode structures 
and quality factors, we shall see in the next section that 
the current results have promising implications. 



VIII. COMPARING MICROSPHERES AND 
FABRY-PEROT CAVITIES 

Fig. [Ii] offers a comparison of the state of the art for 
Fabry-Perot and microsphere cavities for cavity QED, as 
well as projections of likely limits for each. It is interest- 
ing to note that in our projections for the limiting cases 
of each, microspheres allow for a significant improvement 
in the critical atom number, A^o, relative to Fabry-Perot 
cavities. On the other hand, a principal advantage of 
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Fabry-Perot cavities relative to microspheres would seem 
to be significant improvements in the saturation photon 
number, n . The specific specific task at hand would 
then dictate which technology to apply. 

As shown in Fig. |l4|, there has already been some 
progress in coupling atoms to the external fields of a 
microsphere ]3l| . The sphere employed for the work 
of Ref. ]3l|| had a radius of a w 60/im, and quality 
factor Q < 5 x 10 7 , corresponding to a mode volume 
of Vp w 3.7 x 10 3 /im 3 , coupling coefhcient g a /(2n) r* 
24MHz, saturation photon number no = 5.54 x 10~ 3 , 
and critical atom number No = 2.99 x 10~ 2 . If in- 
stead this experiment were to be implemented with a 
smaller sphere with 10/mi radius and with quality factor 
Q ~ 0.8 x 10 7 such as we have manufactured and de- 



scribed in Section VI, the following parameters would 
be achieved: a mode volume of Vp « 1.4 x 10 2 /im 3 , 
coupling coefficient g a /(2ii) w 233MHz, saturation pho- 
ton number uq w 6.27 x 10~ 5 , and critical atom num- 
ber iVo ~ 2.11 x 1CP 3 . Therefore, we see that currently 
achievable quality factors in spheres of radius 10 am al- 
ready would allow for impressive results in cavity QED 
with single atoms. 

By comparison, the state of the art for Fabry-Perot 
cavities has already achieved the following results for the 
TEMnn modes 
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a cavity finesse of T — 4.8 x 10 5 , 
a mode volume of V m « 1.69 x 10 3 ^m 3 , coupling co- 
efficient <7o/(27r) r> 110MHz, saturation photon num- 
ber no 2.82 x 10~ 4 , and critical atom number No ~ 
6.13 x 10~ 3 . If one then looks at possible limits of 
Fabry-Perot technology for cavity QED as analyzed in 
Ref. ^9|, the following may be possible; a cavity of 
length Ao/2 with a cavity finesse of T — 7.8 x 10 6 yields 
coupling coefficient go/(2w) w 770MHz, saturation pho- 
ton number no ~ 5.7 x 10~ 6 , and critical atom number 
N w 1.9 x 10~ 4 . 

It is encouraging that the currently achievable results 
for small sphere manufacture would already allow the 
WGMs to compete favorably with the current state of 
the art in Fabry-Perot cavity QED. However, if one were 
able to manufacture and couple to spheres at the optimal 
size a w 7.83/mi with a Q ~ 9.76 x 10 s (the theoreti- 
cal maximum predicted from the analysis of Section IV), 
the following results could be achieved: a mode volume 
of Vp ~ 90/im 3 , coupling coefficient g a /(2n) « 318MHz, 
saturation photon number no ~ 3.36 x 10" 5 , and criti- 
cal atom number Nq « 9.28 x 10~ 6 . This would repre- 
sent a significant improvement over the current Fabry- 



Perot technology and be competitive with the likely lim- 
its of Fabry-Perot technology. However, even short of 
achieving this stated maximum Q for the WGMs, im- 
pressive results can already be attained. With a qual- 
ity factor Q ~ 0.8 x 10 7 at the optimal sphere ra- 
dius a fa 7.83/xm, one would obtain these same results 
(i.e., g a /(2Tr) rj 318MHz and saturation photon number 
no ri 3.36 x 10~ 5 ), except that the critical atom num- 
ber, No, would increase to No ~ 1.13 x 10~ 3 . This is 
still an impressive gain over the current capabilities of 
Fabry-Perot cavities for the saturation photon number, 
with room for improvement in the critical atom number. 

Overall, we thus find that the technologies of mi- 
crospheres and Fabry-Perot resonators each have their 
advantages and disadvantages. However, there is one 
notable advantage of microspheres; they can be made 
cheaply and relatively simply given sufficient training 
and skill. By contrast, the Fabry-Perot cavities consid- 
ered here require specialized coating runs with expen- 
sive equipment and considerable expertise, which is to 
be found at only a few locations worldwide. This alone 
makes microspheres an attractive alternative to Fabry- 
Perot cavities for cavity QED. Another unique advan- 
tage of the WGMs is the ability to control the cavity 
decay rate, k, by controlling the coupling efficiency into 
and out of the microsphere (e.g., by adjusting the dis- 
tance between a coupling prism and the microsphere 
p2[.) Furthermore, as one moves to the limit of small 
cavities, the open geometry of microspheres offers a con- 
siderable advantage when compared to the geometry of 
Fabry-Perot cavities. Such possibilities combined with 
our projected values of the critical parameters, (no, No), 
shown in Fig. [lj point to the competitiveness of micro- 
spheres with current and future Fabry-Perot technology 
and demonstrate their potential as a powerful tool for 
cavity QED in the regime of strong coupling. 



Acknowledgments 

We thank K. Birnbaum, S. J. van Enk, C. Hood, 
V. Ilchenko, A. Kuzmich, R. Legere, P. Lodahl, T. Lynn, 
H. Mabuchi, J. McKeever, T. Northup, D. Vernooy, and 
J. Ye for helpful discussions. This work was supported by 
the National Science Foundation, by the Office of Naval 
Research, and by the Caltech MURI on Quantum Net- 
works administered by the Office of Army Research. 



[1] V. B. Braginsky and V. S. Ilchenko, Sov. Phys. Dokl. 32, 
307 (1987). 

[2] M. L. Gorodetsky, A. A. Savchenkov, and V. S. Ilchenko, 

Opt. Lett. 21, 453 (1995). 
[3] D. W. Vernooy, V. S. Ilchenko, H. Mabuchi, E. W. Streed, 

and H. J. Kimble, Opt. Lett. 23, 247 (1998). 
[4] D. W. Vernooy and H. J. Kimble, Phys. Rev. A 55, 1239 



(1997). 

[5] H. Mabuchi and H. J. Kimble, Opt. Lett. 19, 749 (1994). 
[6] M. L. Gorodetsky and V. S. Ilchenko, Opt. Commun. 

113, 133 (1994). 
[7] F. Treussart, J. Hare, L. Collot, V. Lefevre-Seguin, D. S. 

Weiss, V. Sandoghdar, J. M. Raimond, and S. Haroche, 

Opt. Lett. 19, 1651 (1994). 



12 



[8] J. C. Knight, N. Dubreuil, V. Sandoghar, J. Hare, 

V. Lefevre-Seguin, J. M. Raimond, and S. Haroche, Opt. 

Lett. 20, 1515 (1995). 
[9] J. C. Knight, N. Dubreuil, V. Sandoghar, J. Hare, 

V. Lefevre-Seguin, J. M. Raimond, and S. Haroche, Opt. 

Lett. 21, 698 (1996). 
[10] W. Jhe and J. W. Kim, Phys. Rev. A 51, 1150 (1995). 
[11] H. Chew, J. Chem. Phys. 87, 135 (1987). 
[12] V. V. Klimov and V. S. Letokhov, Opt. Cummun. 122, 

155 (1996). 

[13] V. V. Klimov, M. Ducloy, and V. S. Letokhov, J. Mod. 
Opt. 43, 549 (1996). 

[14] W. Klitzing, R. Long, V. S. Ilchenko, J. Hare, and 
V. Lefvre-Seguin, New J. of Phys. 3, 14 (2001). 

[15] M. M. Mazumder, D. Q. Chowdhury, S. C. Hill, and R. K. 
Chang, Optical Processes in Microcavities, vol. 3 (World 
Scientific, Singapore, 1996). 

[16] A. J. Campillo, J. D. Eversole, and H.-B. Lin, Optical 
Processes in Microcavities, vol. 3 (World Scientific, Sin- 
gapore, 1996). 

[17] R. E. Slusher and U. Mohideen, Optical Processes in Mi- 
crocavities, vol. 3 (World Scientific, Singapore, 1996). 

[18] A. B. Matsko, S. P. Vyatchanin, H. J. Kimble, and 
H. Mabuchi, Phys. Lett. A 192, 175 (1994). 

[19] D. J. Norris, M. Kuwata-Gonokami, and W. E. Moerner, 
Appl. Phys. Lett. 71, 297 (1997). 

[20] X. Fan, A. Doran, and H. Wang, Appl. Phys. Lett. 73, 
3190 (1998). 

[21] H. J. Kimble, Cavity Quantum Electrodynamics (Aca- 



demic Press, San Diego, 1994). 

[22] H. J. Kimble, Physica Scripta T76, 127 (1998). 

[23] J. A. Stratton, Electromagnetic Theory (McGraw Hill, 
New York, 1997). 

[24] Because the WGMs are the modes of an open resonator, 
the mode volume, Vp, diverges as rg — > oo. However, this 
divergence is logarithmic, and Vp is quite insensitive to 
the choice of rQ for a large range of values. For example, 
at the optimal sphere size discussed in Section N% where 
I = 76, Vp varies by less than 1% for 10 < ^(Tq - a) < 
10 3 . In this case, we take (rq — a) ~ 100/im. 

[25] M. Bass, ed., Handbook of Optics (McGraw Hill, New 
York, 1995), vol. II, chap. 33, p. 63, 2nd ed. 

[26] C. Lin, Handbook of Microwave and Optical Components 
(Wiley, New York, 1991). 

[27] V. V. Datsyuk, Appl. Phys. B 54, 184 (1992). 

[28] H.-B. Lin, J. D. Eversole, and A. J. Campillo, Optics 
Comm. 77, 407 (1990). 

[29] C. J. Hood, J. Ye, and H. J. Kimble, Phys. Rev. A 64, 
033804 (2001). 

[30] C. J. Hood, T. W. Lynn, A. C. Doherty, A. S. Parkins, 
and H. J. Kimble, Science 287, 1457 (2000). 

[31] D. W. Vernooy, A. Furusawa, N. P. Georgiades, V. S. 
Ilchenko, and H. J. Kimble, Phys. Rev. A 57, R2293 
(1998). 

[32] M. L. Gorodetsky and V. S. Ilchenko, J. Opt. Soc. Am. 
B 16, 147 (1999). 



